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Abstract 

Physical collisions and close approaches between stars play an important role in the formation of 
exotic stellar systems. Standard theories suggest that collisions are rare, occurring only via random 
encounters between stars in dense clusters. We present a different formation pathway, the triple evolu- 
tion dynamical instability (TEDI) , in which mass loss in an evolving triple star system causes orbital 
instability. The subsequent chaotic orbital evolution of the stars triggers close encounters, collisions, 
exchanges between the stellar components, and the dynamical formation of eccentric compact binaries 
(including Sirius like binaries). We demonstrate that the rate of stellar collisions due to the TEDI 
is approximately 10""' yr~^ per Milky- Way Galaxy, which is nearly 30 times higher than the total 
collision rate due to random encounters in the Galactic globular clusters. Moreover, we find that the 
dominant type of stellar collisions is qualitatively different; most collisions involve asymptotic giant 
branch stars, rather than main sequence, or slightly evolved stars, which dominate collisions in glob- 
ular clusters. The TEDI mechanism should lead us to revise our understanding of collisions and the 
formation of compact, eccentric binaries in the field. 



1. INTRODUCTION 

The predominant cause of stellar collisions is thought 
to be random encounters between stars in dense clus- 
ters. These collisions, where stars pass within a single 
stellar radius of each other, may result in a wide range 
of outcomes, most notabl y the fo rmation of close and ec- 
centric compact binaries (|Sharalll999l 120021 : lDaviesll2002l : 
[ivanova 2010, and references therein). Because the colli- 
sion rate per star, F, due to random encounters depends 
directly on the local number density of stars F oc n, such 
interactions are common only in the densest stellar en- 
vironments such as globular clusters (GCs). Outside of 
these clusters, the densities are lower by a factor as great 
as 10 : random collisions effectively do not occur in the 
field. In this paper, we describe a dynamical instability 
that occurs in evolved systems containing three or more 
starfQ. This mechanism, which we term the triple evo- 
lution dynamical instability (hereinafter TEDI) , leads to 
close encounters and collisions between stars even in low 
density stellar environments. 

The TEDI occurs due to an interplay between stellar 
mass loss and orbital instability. Consider a hierarchical 
triple system composed of an inner binary coupled to an 
outer binary, where the outer binary consists of a tertiary 
companion in orbit about the center of mass of the in- 
ner pair. When the orbits of the inner and outer pair are 
well separated, the system can remain stable indefinitely. 
Stellar evolution can disrupt this architecture. When the 
more massive component of the inner binary evolves, it 
begins to shed mass. As a result, the orbits expand in 
proportion to the ratio between the initial and final mass 
in th e enclosed system (Hadjidc metriou 1963: EggletonI 
120061) . Since the relative mass loss in the inner binary 
is greater than that in the outer binary, the inner orbit 

^ We focus here on triple systems, but the scenario is also appli- 
cable to higher multiplicity systems. 



encroaches on the outer orbit. This relative orbital ex- 
pansion can trigger a dynamical instability, which results 
in the dissolution of the system (see Fig. [1]). 

The subsequent evolution of such unstable systems is 
very similar to that of binary-single star encounters in 
globular clusters, and often leads to close encounters be- 
tween any of the three components. The TEDI provides 
a unique pathway to physical collisions during these close 
approaches: as stars evolve off of the main sequence, they 
not only lose mass, but also expand radially by a factor of 
a few hundred on the Asymptotic Giant Branch (AGB). 
The simultaneity of mass loss (and thus instability) with 
stellar bloating greatly increases the cross section for col- 
lisions. For this reason, as we now demonstrate, the 
TEDI is the dominant route to stellar collisions in the 
Universe, producing nearly 30 times more collisions than 
random gravitational encounters in GCs. 

The significance of mass-loss induced instability in 
triple systems has been proposed in s everal studies 
(iKiseleva et al.l[T99l llben fc Tutuko-^11999 : Perets 201(| 
iFreire et all 120111 : iPortegies Zwart et al.l i2011i) . How- 
ever, the frequency of this phenomena, and its impor- 
tance for stellar collisions and the consequences thereof, 
have not been explored. 

Our primary goal in this paper is to estimate the rate 
of collisions per galaxy induced by the TEDI mechanism, 
Fcoi. This calculation requires the estimation of the fol- 
lowing quantities: 

• Ni, - The number of stars in the Galaxy 



/, 



evolve 



The observed fraction of evolved stars 



• /triple - The fraction of stars in hierarchical triple 
systems 

• /des - The fraction of hierarchical systems that be- 
come unstable as the stars evolve 



• /col - The fraction of unstable triple systems that 
end with a collision between two of the stars 

Thus Fcol — N^, X /evolve X /triple X /des X /col- 

Observations provide estimates of the first two vari- 
ables, and less certain estimates for the third (triple frac- 
tion). Attaining estimates for /des and /coi occupies the 
remainder of this paper, as well as addressing the un- 
certainties in the triple fraction, /triple- The latter is 
done with both observational estimates and theoretical 
population synthesis models to provide independent es- 
timates. The main difficulties in doing so are due to 
the limited data existing on triple systems and the dis- 
tribution of their properties, and the large phase space 
of possible evolutionary routes. Given these caveats, we 
can only hope to provide an order of magnitude estimate 
of the importance of the TEDI. Though we provide the 
formal statistical uncertainties for the resulting fractions 
we obtain, the systematic uncertainties due to our limited 
knowledge of the distribution of triple system properties 
are likely larger than the statistical ones. 

The outline of the paper is as follows. In Section [2] we 
quantify the criterion for the TEDI. We then determine 
the fraction of systems that undergo the TEDI, /des, in 
Section 131 followed by an estimate of the frequency, /coi, 
of collisions in the destabilized systems O. In Section 
[5] we use these results to provide an estimate of the total 
collision rate in the galaxy. We describe the various con- 
sequences of the TEDI in Section IH and briefly point out 
the influence of secular dynamics on triple stellar evolu- 
tion (!J5]). We close the paper with the discussion and 
summary in Section [3] 

2. THE TRIPLE EVOLUTION DYNAMICAL INSTABILITY 

(TEDI) 

A triple system is dynamically unstable when Q < Qst , 
where Q = a putfl - eout)/ain ai id Qst is (we adopt 
iValtonen et~all 2008; but see also iMardlind 2008 for a 
different definition) : 

Qst = 3(l+m3/Afi2)^/3(i_g^^^)-i/6 Q ^_ 1 cosi 



cos^ i 



(1) 
where ain and Oout are the semi- major axes (SMAs) of the 
inner and outer binaries in the triple system, Cout is the 
outer binary eccentricity, and TO3 and M12 = toi -I- 1712, 
are the masses of the outer component and the inner 
binary, respectively. The mutual inclination between the 
inner and outer binaries is i. 

As the primary (with mass mi) in the inner binary 
evolves off the main sequence (MS), it loses mass. We 
consider only systems where the stellar- wind driven mass 
loss occurs adiabatically, i.e. on a timescale that is long 
compared to the outer orbital period (prompt mass loss 
could occur in supernovae and also lead to instability) 0. 

^ The triple supernova instability scenario is limited only to high 
mass systems. Moreover, in the majority of cases, supernova ex- 
plosions lead to the prompt unbinding of the system (especially if 
they are accompanied by a natal kick) , making longer term chaotic 
evolution, and hence physical collisions, less likely than in the case 
of adiabatic mass loss. The contribution of this process to the over 
all collision rate in the Galaxy is therefore likely to be small, and 
not affect our main results. Nevertheless, this process, and the 
effects of natal kicks, could be especially important for collisions 
with neutron stars and black holes, and merits further study. 



We further assume that mass loss is isotropic, and that 
mass transfer is negligible. The expansion of the inner 
orbit due to mass loss will be larger than the expansion 
of the outer orbit because the relative mass loss in the 
inner system is larger. Thus the SMAs evolve as: 



mi ■ 



- iTi2 TOi,i + m2 + ma 



flout. f 



ain,f 

ain,i mi,f + m2 ' mij + 777,2+ m^ aout,i 



(2) 



()Hadiidemetrioulll963f ) where the subscripts correspond 
to the system components and the evolutionary stage of 
the system [i and / for the initial and the final, post-mass 
loss system, respectively). The expansion of each of the 
SMAs in direct proportion to the mass loss derives from 
the conservation of the quantity GMa{\ — e^), where M 
is the total system mass. Note that slow, isotropic mass 
loss is required for e to remain constant. 

The ratio between the SMAs of the outer and inner 
binaries therefore increases, by a factor of 



1/3 



flin.f/a, 



out,f 



/ 777, 



l,i 



7772 



7771 J +7r72 +?773 

ain,i/aout,i \mi,i^m2) V"7i,i + 7772+7773 



(3) 



The TEDI occurs when this increase drives a stable 
triple system across the threshold defined in Eq. ([T]). 
The subsequent dynamical evolution of an unstable sys- 
tem is chaotic. The system may evolve through var- 
ious configurations, including exchanges between the 
stars in the inner and outer binaries, eccentricity ex- 
citation, and collisions (see Fig. [T] for an example). 
The final outcome is typically the ejection of one star 
and/or a physical collision. Because this three body 
problen i is genuinely chaotic (typical Lyapunov exp onent 
of 1/2: llvanov fc Cherninl Il99ll: iHeinamaki et all 119991: 
IValtonen fc Karttunenll2006l and references therein), the 
outcome is strongly dependent on small variations in the 
initial conditions. 
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Fig. 1. — Three-body simulations of the dynamical evolution of 
triple systems undergoing adiabatic mass loss. The separations (r) 
between each of the three components in each system are plotted as 
a function of time. Also shown are the mass (dashed-dotted) and 
the stellar radius (dashed) of the primary, M\ . Thick lines show the 
pairs including the primary. Top (la): A collision between a MS 
star and an AGB star (the subsequent evolution, assuming point 
masses and no collision is also shown). Middle (lb): the formation 
of an eccentric close WD-MS binary, similar to the Sirius system. 
Bottom (Ic): the WD primary is ejected after several exchanges, 
leaving behind an eccentric MS-MS binary. In all cases the sys- 
tems started with masses of Mi = SM©, M2 = 2.0 M©, M3 = 
2.25 iW0, Oin = 15AU,aout =53AU, eout =0.1, i = 0.1. 



3. ESTIMATE OF THE FRACTION OF SYSTEMS 
DESTABILIZED DUE TO THE TEDI 

To determine whether a given system is unstable, we si- 
muhaneously fohow its orbital and stellar evolution, and 
check whether or not it violates the stability criterion at 



any point during its evolution. We estimate the fraction 
of destabilized systems in two ways: (1) we determin e the 
fraction of observed triple systems ()Tokovininll2008D with 
known orbital parameters that become unstable, and (2) 
we calculate the unstable fraction of a synthetic popu- 
lation of triples modeled on binary statistics to account 
for potential biases in the observed sample. 

We first describe the procedure to identify systems that 
should become unstable, and then calculate the instabil- 
ity rates for the observed and the synthetic samples. 

3.1. Identifying Unstable Systems 

In order to determine whether a given triple system is 
susceptible to the TEDI during its evolution, we make 
use of single and binary population synthesis stellar 
evolution codes (SSE an d BSE, described in detail in 
iHurlev et al.l 120001 120021) , to which we couple an ana- 
lytic calculation of the orbital evolution of the wider third 
component using Equation ([2|). 

The specific procedure is as follows. 

1. The inner binary in the triple is evolved using the 
BSE code up to a Hubble time (13.7 Gyrs), assum- 
ing Solar metallicity. 

2. The third companion is evolved using the SSE code 
up to a Hubble time (13.7 Gyrs), assuming Solar 
metallicity. 

3. The mass-dependent orbital evolution of the inner 
binary and outer single star are combined to form 
a time series evolution of all three stars. The inner 
orbit of the binary is determined by the BSE code. 
The outer orbit is updated at every output time 
step, ti of the BSE using Eq. (2), where the initial 
and final time steps are the preceding and current 
time steps 

O'ovAyti) _ 'T^triplo(ii-l) 



flout (^i-l) TTnf[i-p\c{ti) 

We require that mass transfer between the inner 
binary and the third companion is always negligi- 
ble, so that we may apply the adiabatic formula 
above. If the outer binary separation is smaller 
than 15 AU, we assume non-negligible mass trans- 
fer is possible either following a common-envelope 
phase for the inner-binary, or stellar evolution of 
the tertiary. We therefore consider all triples with 
outer separations smaller than 15 AU to be stable 
to the TEDI. 

4. We calculate the stability coefficient, Q at every 
step in the evolution, and using Eq. (1) find all 
systems that become unstable. Note that if the 
inner binary merges, or any of the stars explode as 
supernovae prior to destabilization, the system is 
considered stable (in the context of adiabatic mass 
loss). All other destabilized triples comprise our 
sample of unstable triples whose size is denoted by 

iVdos. 

5. For all destabilized systems we find the radius, 
i?dcs, of the mass losing evolving star at the point 
of destabilization, and note the types and radii of 
the other components. 



6. The fraction of destabilized systems is defined as 

/dos — -^dcs /triple ■ 

3.2. Observed Triples 

The observed triple sample (|Tokovininll2008[ ) provides 
estimated masses and orbital periods for each star, but 
does not provide their mutual inclinations and eccentric- 
ities. We randomly assign each system relative inclina- 
tion s and eccentricities drawn fro m realistic distributions 
(see lFabrvckv fc Tremainell2007t ). 

We employ a distribution of eccentricities that is 
function of period, in lin e with observed binaries 
(jDuauennov fc Mavoi1ll99H) . For periods shorter than 
1000 days, the eccentricity is chosen from a Rayleigh dis- 
tribution [dp oc eexp(-/3e^)de] with (e^) = /3~^/^ = 
0.33. For periods longer than 1000 days, the eccen- 
tricity is chosen from an Ambartsumian distribution 
{dp = 2ede), which corresponds to a uniform distribution 
on the energy surface in phase space. The inclinations 
are chosen from a random distribution that is uniform in 
cos ?. We repeat the random sampling of these properties 
30 times for each system. 

We only evolve triples that are initially stable, since 
unstable systems are unlikely to be observable due to 
the instability timcscale. If our random assignment of e 
and i produce an unstable system, we draw new values 
from the distributions. 

We follow the procedure discussed in Section [XT] to de- 
termine whether a given system becomes unstable. We 
find that the fraction of all observed triple systems with 
primary mass mi > 1 Mq destabilized within a Hub- 
ble time is /des.o = 0.035 ± 0.006, where the error bars 
correspond to the 1 a statistical uncertainty. 

In Fig. [21 we illustrate the sensitivity of the results to 
the eccentricity distribution: we calculate the fraction of 
unstable systems for a complete grid of eccentricities. For 
this sample we conservatively assume retrograde orbits, 
i = 180°, which are the most stable configurations (see 
Eq. p. 

3.3. Synthetic Sample of Triples 

We complement our use of the existing observed triple 
sample, which may introduce unknown biases, with a 
synthetic samples of triples. The method and the as- 
sumptions used to produc e the syntheti c samples arc 
based on a method used by (iFabrvckv fc T rcmainc 2007) 
for low mass triples and (jPeretsI 120091) for high mass 
triples (with primary mass larger than 3 Mq). 

We assume that the orbital and stellar characteristics 
of triple systems are equivalent to two (uncorrelated) bi- 
nary systems, with the inner binary acting like a single 
mass when choosin g the third component. The study by 
IKraus et al.l (j2011|) suggests that this assumption is in 
good agreement with observations (though this was only 
checked for a small sample of triples). We choose initial 
orbital periods such that both the inner and outer bina- 
ries follow the observed distribution. The inner and outer 
periods are chosen from a distribution that is Gaussian 
for lo w mass binaries with M < 3Mq (jRaghavan et al.l 
120101) and log flat for massive binaries corresponding to 
/(r) oc 1/r i.e., Opik's law (IKobulnickv fc Frveii[2007i) FI 
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Fig. 2. — The fraction of destabilized systems, /dcsi is shown 
for a grid of inner and outer eccentricities applied to the observed 
sample of triple systems (smoothed and interpolated; original bin 
size is Ae = 0.05). 

The mass ratio of the inner binary, q — mi/m2, is cho- 
sen from a Gaussian distribution (mean of 0.6 and dis- 
persion of 0.1) for low mass primaries and a power law for 
massive stars (/(<?) oc q^°*);. The mass of the tertiary 
companion, 7713, is chosen such that q = m^/{mi + m-2) 
follows the same mass ratio distribution. Thus the mass 
of the third star is correlated with the combined mass of 
the inner binary. 

The inner and outer SMAs are computed from the 
masses and periods assuming non-interacting Keplerian 
orbits. The inclination and eccentricity distributions are 
chosen in the same manner as the observed sample. We 
produce a sample of triples with primary masses in the 
range 1 < rrip < 20, with 1000 systems per primary mass 
bin (with bin size AM — 0.3 Mq, for a total number of 
6 X lO'* systems). Note that since the triple outer compo- 
nent is randomly chosen from a distribution with respect 
to the inner binary mass, the term primary mass refers to 
the inner binary primary mass (which is not necessarily 
the most massive component in the triple). 

In Fig. [3^ we show the fraction of destabilized triples 
from this sample as a function of primary mass. Fig. 
[3)d shows the typical radius of the mass losing star at 
the time of destabilization (i?(jes)- The change in the as- 
sumed binary properties between the low mass and the 
high mass stars introduces a sharp change in the typ- 
ical i?dcs, since higher mass stars have more compact 
triple configurations on average, and can be destabilized 
at an earlier phase of their evolution, at which point the 
evolved star has a smaller radius. 

Using the fraction of destabilized triples as a function 
of primary mass, /dcs,s("T')j shown in Fig [31 we can es- 
timate the total fraction of destabilized triples. The 
statistical uncertainty for these fractions is small (1-2 
%; not shown). The fraction of observed stars in triple 
systems ranges from 11% for Solar like (F, G, K) stars 



We assume A stars share a similar distribution as that of F/G/K stars in the l|Haghavan et_aLH2010|) sample. 
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Fig. 3. — Top: The fraction of destabilized triples as a function 
of the inner binary primary component mass. The down turn at 
about 8 M0 is due to supernova explosions that occur before the 
system can be destabilized due to mass loss. Bottom: The median 
radius of the mass losing stars at the onset of instability, as a 
function of mass. Separate lines indicate the low mass (rra < 3 
Mq) and the high mass stars (m, > 3 Mq), for which different 
triple distributions are assumed (see text). 



(jRaghavan et al.ll20ia lRastegaevll2010D to possibly 50% 
for more massive B stars (|Evansll201ll) 7^We use a Salpeter 
initial mass function, g{m) = dn/dm oc m~^'^^, and as- 
sume the triple fraction of massive s tars is about 5 times 
large r than that of lower mass stars (|Kobulnickv fc Frveii 
I2007t lEvansll201lD . Given these distributions, we esti- 
mate the total fraction of destabilized triple systems in 
our synthetic sample to be: 



/dcs,s = /i /dcs,s(rn)^rf"^ 



r-30 



+5 X /g"" /des.sM^dm = 0.010 + 0.043 = 0.053.(4) 

The overall fraction of destabilized triples in our syn- 
thetic sample is higher than but consistent with the 
results for the observed sample (0.053, compared with 
0.035) (see also Fig. [2]). Our estimates suggest that 
a few percent of all evolved triple systems undergo the 
TEDI. 

We now turn to direct integrations of three-body sys- 
tems to both validate the analytic instability criterion 



and to determine the fraction of unstable systems that 
produce a collision, fcoi- 

3.4. Three-body Simulations 

Because TEDI systems begin in stable configurations 
and adiabatically evolve into unstable configurations, 
they may not be directly compar able to unstable triples 
explo red in previous studies (jValtonen fc KarttunenI 
|2006[ ). Therefore we conduct a series of three-body simu- 
lations including mass loss both to check that our desta- 
bilized systems are truly unstable, and to estimate the 
frequency of collisions based on the statistics of close ap- 
proaches during the chaotic evolutionary phase. In the 
appendix, we provide a semi-analytic estimate of the col- 
lision rate. 

The integrations are done using a modified ve r sion o f 
the Hermite integrator described in iHut et ahl ()1995D . 
The integrator uses a variable time step constrained to 
be 10~^ of the minimum collision time of any two stars. 
Reducing the time-step by a factor of 10 does not change 
the outcome. In runs without mass loss, energy is con- 
served to about one part in 10^" over the 10 Myr run 
time. We only consider point masses, and include no 
tidal effects, or mass transfer. We use a constant mass 
loss rate, determined by the total mass loss in the system 
as found in the stellar evolution calculation, divided by a 
constant mass loss timescale, rioss = 0.5 Myr (the typical 
lifetime of the highly evolved stars, during which most of 
the mass loss occurs). After the mass loss epoch, rioss, 
we continue the simulation, at constant mass up to 10 
Myrs to check for longer term stability. In cases where 
our stellar evolution calculation shows that instability 
occurred only after both the primary and secondary in 
the inner binary evolved, we simulated the system begin- 
ning only at the second stage of mass loss (i.e. from the 
secondary) , after the primary became a white dwarf. 

We simulated 100 realizations of each observed triple 
system found to be unstable in SectionlH each realization 
differed only in the initial orbital phase of the stars. Since 
we do not account for mass transfer, we simulated only 
systems in which the inner binary separation was larger 
than 15 AU (these comprise nearly half of the observed 
sample). We ran 300 different system configurations, re- 
sulting in a total of 3 x 10'* runs. 

To identify systems that become unstable, we require 
that either one star is ejected from the system, or an ex- 
change occurs between the tertiary and one of the inner 
binary components. We find that approximately 0.33 of 
the systems found to be unstable in our simplified cal- 
culations become unstable in our simulations over a 0.5 
Myr run. Although more triples are likely to show ex- 
changes and escapes over longer time scales (for example 
we find approximately 55% were destabilized within 10 
Myrs), we adopt the short timescale simulation result, 
since our primary concern is with collisions, which are 
much more likely if they occur while the primary is on 
the AGB. 

We take the most conservative estimate possible by 
assuming that the results from the three body inte- 
grations are more representative of the entire sample. 
We thus reduce the fraction of unstable observed sys- 
tems found in Section [3] (/dcs,o = 0.035) by a factor of 
0.33; we adopt as the fraction of destabilized systems 
/des = 0.035 X 0.33 ~ 0.012. 



4. STELLAR COLLISIONS DURING THE TEDI 
(CALCULATING Fcol) 

Having calculated the expected fraction of triple sys- 
tems that are unstable to the TEDI, we now estimate the 
fraction of such systems that might suffer physical colli- 
sions. We show in the appendix that applying the results 
of previous simulations for both single star - binary en- 
counters and unstable triple systems gives the fraction of 
unstable systems that undergo a collision as 0.49 ± 0.03. 
The analytically derived collision rate for unstable triples 
is significantly higher than that for random encounters 
because an unstable triple typically undergoes multiple 
close approaches before the system dissolves, and because 
the stellar radius at the time of collision is so large. Be- 
cause the analytic results are based on initially unstable 
systems covering a limited part of the the three body 
problem phase-space, we determine a more conservative 
collision rate based on our three body integrations (see 
below). 

4.1. Collisions in three-body simulations 

We define collisions in our three body integrations as 
those cases where the closest approach between the pri- 
mary and another star is smaller than the typical stellar 
radius of the evolving star, i?des where i?dcs is defined 
to be the stellar radius at the point of instability, when 
the system first satisfies Q < Qst- We only classify close 
approaches as collisions if they occur within the 0.5 Myr 
mass loss timescale when the star is likely bloated. Figure 
m shows the cumulative distribution of closest approaches 
amongst the unstable triples that we simulated. The dis- 
tribution of radii of the evolved stars that go unstable 
(i?dcs defined before) is also shown. Most stars involved 
in the collisions are highly evolved AGB stars, with radii 
in the range of tens to hundreds of AU. 

We find the collision rate between the evolving star 
and either of the other stellar components to be /coi — 
0.09 ± 0.03 amongst the unstable systems, where the er- 
ror bars are the 1 a statistical uncertainty. Note that 
this fraction is smaller than that obtained from the sim- 
plified semi- analytic approximation discussed in the ap- 
pendix. However, detailed analysis of the distribution of 
closest approaches between stars during the triple evo- 
lution, suggests that the simple linear approximation in 
Eq. ()A1|) is inaccurate; we find that the closest approach 
statistics depend more complexly on the various orbital 
properties of the system. We therefore adopt the colli- 
sion fraction found in our simulations for our subsequent 
calculations. 

The sample chosen for the three body integrations ex- 
cluded close inner binaries (with SMAs smaller than 15 
AU), because we did not account for mass transfer. Phys- 
ical collisions are even more likely for these systems due 
to the monotonic dependence of the close approach prob- 
ability on the inner binary separation (seen both in our 
integrations and in Eq. lAll ) Thus the neglected sys- 
tems likely have a higher collision probability. Because 
we adopt the same collision fraction for more compact 
triples, these rates should be considered a lower bound. 
Compact systems comprise about half of our observed 
sample, so the neglect of these systems in our calcula- 
tions would change our results by at most a factor of 
two. 
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Fig. 4. — Stellar radii of evolved stars, and closest approaches 
between stars in destabilized systems. Dashed red lines indicate 
the distribution of closest approaches between the evolved star and 
another member of an evolving triple during the mass loss phase 
(taken to be ti^ss = 0.5 Myr). Solid blue lines show the distribution 
of the stellar radii of the evolved mass losing star at the point of 
instability. Inset shows the cumulative distributions. Results were 
obtained from the three-body simulations of systems that became 
unstable. 

Note that we cannot correctly follow the full evolution 
of systems in which tidal effects dominate the dynam- 
ics at any point. Nevertheless, for the short timescale 
of the TEDI, tidal effects become important when the 
closest approach is comparable to the stellar radius, and 
can be considered a sub-type of collision. Moreover, be- 
cause most collisions occur when the physical radius is 
comparable to the tidal disruption radius, our neglect of 
tidal effects should not significantly alter our calculated 
collision rates. Determining the configurations of post- 
collision systems is beyond the scope of this work. 

5. GALACTIC COLLISION RATE DUE TO THE TEDI 

In the previous sections we provided estimates for the 
fraction of destabilized systems, /dcs, and the fraction 
of colliding systems, /coi, among these systems. We can 
now make use of these results to estimate the overall rate 
of stellar collisions due to the TEDI. Taking the estimate 
of the collision rate together with the calculated fraction 
of destabilized triples, we find that the collision rate be- 
tween evolved and companion stars/compact objects in 
a Milky Way like galaxy, over the age of the Universe, ^h 
is: 

T^ _ ^i< ^ ,/ evolve X /triple ^ /des ^ Jeol ,-x 

i eol — ; — (,0J 



tH 



= 1.2 X 10" 



TV, 



1.5 X 10" 



/evolve \ / /triple \ / /des \ / /col \ _i 



0.1 



0.1 J V 0.012/ V0.09 



where N^, is the n umbe r of stars in the Galaxy 

the mean stellar 

/evolve is the ob- 

served fraction of evolved stars (the fraction of WDs, 



(jBinnev fc Tremaind l2008t taking 
mass to be 0.3 Mq; iKroupal [20011) . 



iBelczvnski et al.ll2002l can serve as a proxy for their total 
number), /triple is the fraction of triple systems (where we 
conservatively adopt the tr iple fraction of / triniR = 0.1, 
found for F/G/K stars, Rag havan et alJ[201ClL for stars of 
all masses). We adopt /dos (the fraction of triple systems 
that become unstable) and fcoi (the fraction of unstable 
systems in which a collision occurs) from our calculations 
in the previous sections. 

For comparison, the collision rate in the GC M15, 
one of the densest in the galaxy {uc > 10^ M0 pc~^; 
lUmbreit et al.ll2008l) was es timated to be F = 2^1x10"^ 
vr~^ (lUinbreit et all [200I see also refs. iSharal 119991 : 
iLee et al.l 120101 ) The total collision rate in Galactic 
GCs, which is dominated by such clusters (i.e. over 
100 GCs exist in the Galaxy, but the dominant contri- 
bution to coUsions come from the most dense, massive 
clusters), is therefore Fgai = 4 x 10"^ (iVcc/SO) yr^^, 
where Nqc ~ 20 is an estimate of the number of sim- 
ilar GCs in the galaxy as massive and dense as Ml 5 
(jGnedin fc Qstrikeilll997[ ). Collisions due to destabilized 
triples are therefore about 30 (A^gc/20) times more fre- 
quent than collisions due to random encounters in GCs. 
The number of extragalactic GCs is typically propor - 
tional to the host galaxy mass (jBrodie &: Stradenl2006l ). 
suggesting a similar ratio of field to GC collisions in other 
galaxies. 

Note that the collsion rate due to TEDI inside GCs 
scales with the total number of stars in GCs and does not 
depend on their specific density. Taking approximately 
150 GCs exiting in the Galaxy with about 10^ stars in 
each cluster, we find the total collsion rate in GCs due to 
TEDI to be at most 0.1 percent of the total Galactic rate, 
and only about 3 percente of the collsion rate in GC due 
to random encounters. In other words, inside GCs, sel- 
lar collsions are dominated by random encounters rather 
than by triple stellar evolution. 

6. CONSEQUENCES OF THE TEDI 

The TEDI scenario has a wide variety of possible out- 
comes. Many of these are very similar to those resulting 
from binary-single and single-single close encounters be- 
tween stars and/or compact objects that occur in glob- 
ular clusters. However, the TEDI scenario provides a 
pathway to these exotic scenarios previously reserved for 
dense stellar populations. 

In the following we discuss various outcomes of the 
TEDI and highlight their unique aspects, either in terms 
of their qualitative difference from typical encounters in 
clusters, or their specific relations to systems in the field. 
We first briefly overview several potential outcomes, and 
then discuss in more detail a specific example, providing 
a novel formation scenario for Sirius like eccentric WD 
binaries. 

Stellar collisions in the field: As discussed above, 
the TEDI produces a high rate of stellar collisions in the 
field, most of which involve AGB stars; these were pre- 
viously thought to comprise a negligible fraction of all 
collisions. These encounters might be observed as in- 
termediate luminosity optical transients (with total en- 
ergy ranging between lO''^^ — lO''^ ergs; these could be de- 
tectable in current and fu ture optical transient surveys; 
iKulkarni fc Kasliwalll2009[ ). possibly similar to those sug- 
gested to occur due to mergers or thr ough tidal/off 
axis collisions in eccentric binaries (jKashi fc Sokcr ,2010; 



lSmith|[201lD . Beside potentially producing immediate 
energetic transient events, the TEDI can produce a wide 
variety of collision products, including merged stars and 
blue stragglers as well as peculiar binaries. 

Numerical simulations of close encounters with stars 
bloated to 100 Rq show that the resulting collisional 
and/or tidal capt ures lead to the formation of highly 
eccen tric binaries (jBailev fc Davieslll999l : lYamada et al.l 
I2OO8I ). Such collisions might also lead to a common en- 
velope phase and the formation of a close binary. The 
TEDI scenario therefore predicts the existence of a sub- 
population of close, but eccentric, evolved binaries in the 
field. The TEDI can explain the puzzling orbits of WD 
binaries such as the nearby Sirius system (see below), 
as well as highly eccentric Barium binary systems (ex- 
pected to form through mass transfer in a close binary). 
These systems pose a problem for standard binary stellar 
evol ution models, which produce on l y circular close bina- 
ries (iBonacic Marinovic et al.ll2008t iKarakas et al.ll2000l: 
llzzard et al.ll2010D I'but see iDermine et al.l ('20"l2l. and 
references therein, for some alternative suggestions). The 
peculiar close WD-blue straggler binary (KOI- 74) found 
by Kepler, which may hav e formed in a h ighly eccen- 
tric evolved binary system (|Di Stefanoll2011l ). could eas- 
ily be produced by the TEDI. Similarly, TEDI evolution 
can lead to the formation of eccentric WD- WD binaries, 
which could provide unique sources for future gravita- 
tional wave detecti on missions, and we re predicted to 
exist only in GCs (Wi llems et al.l [20071 ). Our findings 
suggest that there could be a few times more such sources 
in the field than estimated to exist in GCs. 

Neutron stars/black hole binaries: Mass loss 
from primaries more massive than 8 M© which undergo 
the TEDI, might form a close eccentric binary before the 
supernova explosion of the evolving massive star. Since 
such binaries have a high er probability of surviving the 
SN explosion ()Hillslll983[ ). the TEDI might aid in the 
retention of neutron stars and black holes in binaries, 
and the later formation of X-ray and pulsar binaries. 
Such triples could later evolve, during mass loss/mass 
transfer from the other triple components. Possible out- 
comes include instability and thus collisions with a neu- 
tron s tar/bl ack hole, or even exchanges. S ee iFreire et al.l 
(|201H) and iPortegies Zwart et al.l (|201lD for the latter 
possibility, but note that these studies did not account for 
the possibility of a physical collision during the chaotic 
evolutionary phase. 

Envelope stripping through collisions: A collision 
between a neutron star and its evolving companion could 
result in the tidal disruption of the latter. The rem- 
nants of this star might form an accretion disk sur- 
rounding the neutron star, which could provide the nec- 
essary mass and angular m omentum to create a recy- 
cled pulsar (|Verbunt et al.l[r987j . If the third compan- 
ion is ejected during the course of the instability, the 
pulsar would appear as an isolated recycled pulsar in 
the field. If the third companion is retained it will 
be observed as a wide and likely eccentric recycled bi- 
nary pulsar in the field. The latter scenario provides 
an additional channel for the forrn a tion of the pulsar 
J1903-H0327 (iChampion et al.l [20081: iFreire et al.l I20TII: 
IPortegies Zwart et al.ll201lD . Another possibility is the 
formation of a recycled pulsar with a stripped compan- 
ion on an eccentric orbit, if the companion's core stays 



intact. 

There are many observed peculiar stars and com- 
pact remnants thought to form through envelope strip- 
ping by a close c ompanion, such as hot sub-dwarf stars 
JHan et alj I2002D . helium white dwarfs (Mars h et aTl 
1995) , type Ib/c supernova progenitors ( Yoon et al.l 
20101 ). However, cases are found where these peculiar ob- 



jects are observed to be single or to have a wide compan- 
ion, configurations that are inconsistent with the close bi- 
nary progenitor scenario. Such cases might be explained 
through a TEDI induced collision. Further studies are 
required in order to quantitatively explore each scenario. 

7. FORMATION OF THE SIRIUS WD-BINARY SYSTEM 

One of the closest WD binaries to Earth is the Sirius 
system composed of a MS star of mass 2 Mp, and a WD 
companion of about a Solar mass (|Liebert et al.l [20051 ). 
The SMA of 20 AU and and high eccentricity of 0.59 
make this binary quite peculiar, as standard stellar evolu- 
tion scenarios suggest that the progenitor of such a close 
binary would have circularized during stellar evolution of 
the system (|Bonacic Marinovic et al.l 120081 : llzzard et al.l 
12010). However, the Sirius system could have formed 
through the T EDI; alter native scenarios are also in devel- 
opment ( ,Dermine et al.l 2012; P. Eggleton, private com- 
munication 2012). 

To assess the likelihood that Sirius formed via this 
mechanism, we ran 6500 three-body simulations of the 
dynamical evolution of initially stable triple systems un- 
dergoing mass loss. For these systems, we follow the ex- 
act mass loss and radius evolution as obtained from the 
SSE code. The primary starts out with M = 5.05 M©, 
and ends as a 0.98 Mq WD. We evolve the systems ana- 
lytically until the primary has reached about 4 Mq and 
begins the rapid mass loss phase. We assume that the 
secondary and tertiary components have undergone no 
evolution. We explore a range of inner binary separa- 
tions from 15 — 30 AU, outer binary separations from 
3.5 — 5.5 times the inner separation, eccentricities from 
— 0.25, and companion masses from 0.6 — 5.5 Mq. One 
component is always constrained to be 2 Mq , the mass of 
Sirius A. The minimum separation of the inner binary is 
chosen to avoid mass transfer. To account for variations 
in the chaotic evolution at late times, we run each set 
of parameters starting from five different initial orbital 
phases. We limit the initial inclination of the systems to 
i = 0.1 radian. 

In Fig. [21 we show examples of final binary configu- 
rations of unstable systems where one of the stars was 
ejected and a Sirius like binary was left behind. Also 
shown is the observed SMA and eccentricity of the Sir- 
ius system. The Sirius binary falls within one of three 
populated regions of final configurations, suggesting that 
Sirius could have formed through the TEDI process. The 
inner cluster of eccentric systems is formed via ejection 
of the wider tertiary companion, while the outer cluster 
of eccentric systems formed via an exchange followed by 
the ejection of the initial inner binary companion to the 
WD. The moderately eccentric systems are composed of 
bound triple systems that underwent only mild eccen- 
tricity excitation. 

We note that if Sirius and other objects formed 
through this mechanism, it might be possible to iden- 
tify its ejected third star, as a very wide companion 




Semi-Major Axis (AU) 

Fig. 5. — The initial and final eccentricities and semi-major axes 
of triple systems evolved using our three-body simulations. Circles 
represent the initial orbital configuration of a given pair of stars 
in the triple. Note that each circle corresponds not only to multi- 
ple realizations of a given system, but also to different masses for 
the third component. Lines connect the initial and final configu- 
rations, which are marked by stars. Each initial configuration can 
correspond to multiple final outcomes due to the chaotic nature of 
the instability. Final configurations similar to the observed Sirius 
system (marked with a red triangle), or even more eccentric are 
a typical outcome for systems with initial SMAs in the range 15- 
30 AU. Also shown are binaries formed via an exchange between 
the inner and outer companions, followed by ejection of the initial 
inner companion. 



(e.g. iJiang fc Tremaiiigl2010( ) or through similar meth- 
ods used to locate stars that originally formed clos e to 
the Sun (jPortegies Zwartil2009l: fBobvlev et al.ll201lD . 

8. TRIPLE STELLAR EVOLUTION COUPLED WITH 
SECULAR KOZAI CYCLES 

Triple systems in which the inner and outer orbits 
have a high mutual inclination may secularly evolve 
through Kozai cycles, in which the inclination and ec- 
centricity of the inner binary can period ically /quasi- 
periodically change by a large amplitude (jKozail [19621 : 
lLidov|[T962l). and even fl ip from prograde to retrograde 
orbits ()Naoz et al.l 120111 : and vice versa) . It was al- 
ready suggested that the coupling of secular Kozai evo- 
lution with dissipative processes such as tidal friction or 
GW emission could play an important role in the cvolu- 
tion of the inner binaries (e.g. [Mazeh & Shaham 1979 



Eggleton fc Kisseleva-Eggleton|[2b06; .Pcrcts fc Fabrvckvl 
20091 : lThompsonll2011l and references therein) . Mass loss 



evolution in triple systems could also couple with this sec- 
ular Kozai process, and change its outcomes. Our simula- 
tions suggest that mass loss does not quench Kozai cyles. 
The high eccentricities that could be induced through the 
Kozai mechanism can therefore lead to close encounters 
between the components of the inner binary during peri- 
center approach prior to destabilization during the AGB 
stage. Kozai evolution could lead to close interactions 
between triple components at a late epoch in the evolu- 
tion, even if the original triple did not interact while 
on the main sequence (see also iShappee fc ThompsonI 
2012). The coupled stellar-dynamical evolution TEDI 
stable, Kozai susceptible, systems is beyond the scope of 
this paper, and will be explored elsewhere. 



9. DISCUSSION AND SUMMARY 

In this paper we have studied the triple evolution dy- 
namical instability (TEDI), in which mass loss in an 
evolving triple system leads to orbital instability. The 
chaotic evolution that ensues gives rise to close encoun- 
ters between the stars and/or exchanges between them. 
Such evolution typically ends in a collision between two 
of the stars or the ejection of one of the components. This 
instability can lead to the dynamical formation of eccen- 
tric compact binaries, and produces a high stellar colli- 
sion rate between stars in the field. The TEDI scenario 
results in a collision rate in the field of approximately 
10~^ yr~^ per Milky- Way Galaxy, about 30 times higher 
than the collision rate in all Galactic GCs. This high 
rate is attributable to two features of the TEDI. First, 
the chaotic orbital evolution results in binary-single star 
like resonant encounters in the field, which provide nu- 
merous opportunities for collisions. Secondly, the cross 
section for physical collisions is much larger than in typ- 
ical random encounters. The mass loss that precipitates 
the instability reaches its maximum at the AGB phase, 
and therefore the instability typically occurs when the 
stellar radius is at its peak, a few hundred Solar radii (as 
seen in Fig. 14]). The cross section for collisions is there- 
fore increased proportionally. In contrast, random colli- 
sions in GCs arc heavily dominated by MS and sl ight! 
evolved stars ([Lombardi et al.ll2006l : lUmbreit et al.ll200 
with radii of 1 — 10 Solar radii. The discrepancy between 
the peak collision phases is due to the 10^ — 10^ year 
duration of the AGB phase. This timescale is a tiny 
fraction of a star's life, and thus short compared to the 
random encounter timescale, even in the densest environ- 
ment. In contrast, the close encounters during the TEDI 
always coincide with the AGB phase. The TEDI model 
therefore requires not only a revision of the rate of stel- 
lar collisions, but also a reconsideration of their nature; 
most stellar collisions involve AGB stars. 

In addition to enhancing the number of stellar col- 
lisions, the TEDI suggests that various types of ex- 
otic binaries collisionally formed inside GCs (X-ray bi- 
naries, cataclysmic variables, single and binary recy- 
cled pulsars, eccentric WD binaries, including eccentric 



ins piraling WD- WD GW sources, and stripped stars; 
F abian et al.l 119751: iPoolev et al.l I20Q3I: iBailev fc Daviesl 



Il999t iLombardi et all 120061 iWillems et all 120071 ) can 
similarly collisionally form through TEDI outside of 
dense clusters; adding an additional evolutionary chan- 
nel for the formation of such binaries in the field. In 
particular, the TEDI provides a novel explanation for 
the origin of WD binaries observed in puzzling configura- 
tions, such as the well known Sirius system. The exotic 
stellar systems produced through these pathways have 
an outsize impact on our understanding of fundamental 
physics, and the evolution of stars in the field. They serve 
as the main sources for gravitational wave emission, and 
as unique probes of the structure of compact objects. 

This study focused on the TEDI scenario for triple 
stellar systems. However, a system comprised of a stel- 
lar binary and a planet on a circumstellar orbit around 
an evolving star can also suffer a similar fate (jPeretsI 



120101 120111 ). The e volution of such sy s tems has recently 
been described by iKratter fc PeretsI (|2012[ ). In addi- 
tion, the TEDI scenario can be extended to higher mul- 
tiplicity systems (note that 1/4 of all multiple systems 
are q uadruples or higher order multiples Raghavan et aLl 
120101 ). Such systems can produce an even wider variety 
of outcomes than triple systems, and in general mimic 
almost any type of outcome that is expected to occur 
following close encounters in GCs. 
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APPENDIX 

SIMPLIFIED ESTIMATE OF THE RATE OF STELLAR COLLISIONS IN THE TEDI SCENARIO 

The frequency of close approaches for initia lly unstable systems has been characterized in pr evious studies (see 
iSaslaw et al.lll97J : iValtonen fc Karttunenl[2006l and references therein; see also lSigurdsson fc Phinn cv 1993 for the re- 
lated issue of binary-single encounters) . Triples starting from unstable configurations undergo simi lar evolution to that 
observed in simulations of low velocity, resonant encounters between a binary and a single star dValtonen fc KarttunenI 
I2006D . The evolution of unstable triples can thus be modeled as a set of random encounters between an incoming 
third body and an inner binary. After any encounter, the incoming star can be ejected from the system, go through 
an exchange with one of the binary components, or collide with one of them. If a collision or ejection does not occur 
on a given passage, the cycle repeats, and the system experiences another encounter (though the components of the 
inner binary may change in the event of an exchange). These encounters cease when a collision or ejection occurs 
(a long-lasting stable triple can only be formed via energy dissipation, e.g. following a collision or a strong tidal 
encounter,). Previou s three body integration s show that, on average, approximately 60 such encounters occur before 
a star is ejected (see IValtonen fc KarttunenI [2006i and references therein), thus the probability for a close approach 
between any two stars in the system is proportionally increased compared with single-single stellar encounters. These 
results provide a simple estimate for the collision rate. During the chaotic ev olution of unstable triple sys tems, the 
probability, P(< rdosc, flout) for a close approach of rdosc ^ flout or smaller is (jValtonen fc Karttunenl[2006[ ): 
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where < L/Lj,jax ^ 1 is the normalized angular momentum of the system {Lmax is the angular momentum for a 
circular orbit of the same energy) . For the typical system that becomes unstable, Eq. (jAll) gives a probability of order 
unity because the stars are often lOO's of solar radii at the onset of instability; however, the time spent at this radius, 
^R.dos, may be shorter than the typical length of the chaotic evolution, tchaotic ~ 60 x PoJjt, where p^"^ is the period 
of the triple outer orbit at the point of destabilization. Thus the probability for a collision is quenched by a factor of 
approximately min(tR4os/60poutJ !)• 

Applying this correction factor, we can calculate the collision fraction for our observed sample. We take Rdcs to be 
the stellar radius of the AGB star at the point of destabilization, Oout to be the outer binary separation, and ^r^os 
to be the AGB lifetime at this radius or above. All these parameters are calculated separately for each system in our 
observed triples sample, at the point of instability, using the models described above (Section [33} ■ 

According to this calculation, the fraction of destabilized triples that suffer a direct collision between the evolving star 
and one of the other triple components is 0.49 ± 0.03 for the observed sample of triples (with the assumed eccentricity 
distribution). Comparison with our numerical results suggests that the analytic calculation significantly overestimates 
collisions for TEDI systems. 
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